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ABSTRACT 



A theory of vector gravitational waves is presented on the basis of an addi- 
tional gravitational field, recovering ideas that go back to H. A. Lorentz, field that 
A. Wheeler called "gravitomagnetism" . The theory is independent of Einstein's 
General Relativity Theory, but is consistent with the Einstein-Lorentz frame- 
work of Special Relativity, so that the resulting equations are considered exact 
from the start, and not approximations valid for velocities small compared with 
the speed of light. A simple model of mass formation and presence of angular 
momentum started with a new type of singularity, "trio-holes" , is exhibited. 



1. Introduction 



The main theoretical purpose of the present work is to exhibit gravitation as a 
particular form of interaction between matter and light (and more generally radiation), 
within the framework of Special Relativity, but in a manner wholly independent of Einstein's 
General Theory. 

The mathematical modeling of gravitational phenomena has not reached yet the state 
of development provided for electromagnetic phenomena by the Maxwell-Lorentz system 
of linear p artial differential equations and its furth er quantization by Dirac, Feynman and 



others (see 



Feymann &: Hibbs 1965 



Feyman 19851 ). After Faraday introduced the notion 



of an electromagnetic field as a physical entity on its own right, the classical theory of 
electromagnetism, in spite of Coulomb's inverse square distance Law of force, could be 
considered free from an action at a distance postulate only thanks to Maxwell's equations. 



In Einstein's own words ( lEinsten 1950l ): "It was Maxwell who fully comprehended 



the significance of the field concept; he made the fundamental discovery that the laws of 
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electrodynamics found their natural expression in the differential equations for the electric 
and magnetic fields. These equations implied the existence of waves, whose properties 
corresponded to those of light as far as they were known at that time.... The new insight 
made it possible to dispense with the hypothesis of action at a distance, at least in the realm 
of electromagnetic phenomena; the intermediary field now appeared as the only carrier 
of electromagnetic interaction between bodies, and the field's behavior was completely 
determined by contiguous processes, expressed by differential equations.... Because one 
cannot dispense with the field concept, it is preferable not to introduce in addition a carrier 
with hypothetical properties." 

Lorentz, did much to extend Maxwell's work into the atomic scale, i ntroducing th e 
Lorentz force and explained the normal Zeem an effect for the first time ( iLorentz 19091 ). 



According to L. De Broglie (jPe Broglie 19391 ) a physicist image of reality needs to change 



when transporting it from the human or astronomical scale onto the much smaller atomic 
scale, where Quantum Laws cannot be avoided so far. In this first work we shall not 
consider the analog of the quantization procedure applied to electrodynamics, considered 
by Dirac and followers, leaving aside the atomic scale. 

There are many reasons to hold on to the group of Lorentz transformations for as 
long as doing this yields fruitful results. Seemingly, Maxwell's equations are here to stay 
at least in the human and astronomical scales, and as Lorentz discovered early, they are 
covariant (for empty space) with respect to Lorentz transformations. While it is true, as 
Einstein observes, that "Maxwell's equations imply the Lorentz group, but the Lorentz 
group does not imply Maxwell's equations", it is also true that "... it is to be expected 
that all equations of p hysics are cov ariant with respect to Lorentz transformations (special 



theory of relativity)" fiEinsten 1950l ). 



Einstein further declares that: "A consistent field theory requires continuity of all 
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elements of the theory, not only in time but also in space, and in all points of space. Hence 
the material particle has no place as a fundamental concept in a field theory. Thus even 
apart from the fact that gravitation is not included. Maxwell's electrodynamics cannot be 
considered a complete theory." 



We doubt that such a statement, dating back to 195 0, would have 



Deen written after 



the advent of Laurent Schwartz's theory of distributions (ISchwartz 19591 ). one of the major 
mathematical breakthroughs of the 20th century, together with Lebesgue's measure and 
integration theory (1901) and Hilbert and Banach's infinite dimensional normed vector 
spaces, among other mathematical advances pertinent to mathematical physics (of course, 
pertinency changes with time). Distribution Theory gave a precise mathematical meaning 
to the familiar Dirac's "delta function 6" (introduced intuitively in the nineteen-thirties) , 
defining it as a continuous linear functional acting upon a topological (non normed) infinite 
dimensional vector space of well behaved test functions, thus allowing for the correct 
handling of partial distributional differentiation of all orders, and also making full sense of 
A(l/r) = —Ait 6 as an equality among distributions ("r" denoting the distance of a variable 
point in £"3 to the origin 0). 

Schwartz's original work was technically difficult, involving the theory of topo- 
logical vector spaces, but afterw' ards, Vladimirov's and Zemanian's expositions (see 



Vladimirov 1984 



Zemanian 19871 ) made the subject accessible to all scientists dealing with 



linear differential equations. 

In the case of non-homogeneous linear partial differential equations (PDE) with 
constant coefficients, taken in a distributional sense, and with a Dirac's delta distribution 
used as a forcing term, the mathematically singular fields obtained by solving them, called 
"fundamental solutions", can be usefully associated with the physical fields generated by 
physical monopoles, as exemplified by the fundamental solution for the Laplace's operator. 
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The distributional derivatives of Dirac's delta adequately symbolize multipoles, the delta 
itself standing for a monopole. Due to the linearity of the distributional differential 
equations, the distributional solution when the forcing term is a multipole is obtained 
from the fundamental solution by distributional differentiation. Besides, the fundamental 
solution can always be obtained by m eans of the distributional Fourier Transform when it 



cannot be obtained by inspection (see IVladimirov 19841 : IZemanian 19871 ). 



In regions of regularity, the distributional solution obtained can be identified with 
an ordinary differentiable function, recovering classical solutions. Just about the only 
thing you cannot do with distributions is to multiply them in a manner that would be 
analogous to the familiar non-linear terms of classical PDE's and fields, except in some 
particular cases not devoid of interest for a physicist, such as those worked out by two of 
the present authors in connection with a uniqueness problem set up by Fritz John for the 
non-linear conservation la w defined by Ut + uu^ = in the domain t > 0, — oo < a; < oo 



( iMiranda fc Pereyra 20071 ). and solved by means of a "moving" delta starting from (0,0), 



with an intensity varying with time. 

Should a linear field theory for gravitation be amenable to a reasonable mathematical 
justification, and also to experimental verification in either the astronomical scale or human 
scale laboratories at the beginning, and also at the atomic scale after some adequate 
quantization procedure had been completed, then by the distributional handling of the 
proposed linear PDE's for adequately defined gravitational fields, it would have been 
established a much simpler way of dealing with known gravitational facts (or allegedly so) 
and also with unexplained features related to angular momentum in the astronomical scale. 
This is the plan we propose to follow in the present work. Besides, the way the theory 
is reasoned, makes of it a more natural step towards a unified treatment of gravity and 
electromagnetism. 
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Einstein sought to unify electromagnetism and gravitation by introducing a curved 
space-time geometrization of the gravitational field first (general theory of relativity), and 
then trying to geometrize electromagnetism. Our point of view will be the opposite, namely, 
we will "electromagnetize" (so to speak) gravitation. 

There have been other theories of gravitation prese nted to the scientific commu - 



nity with varying evidence of experiraental support (see 



AUev fc Yilmaz 1995 



Mashhoon 2003 



Hermann 1978 



AUev 1995 



Salam et al. 19981 ). and Einstein's 1928 work 



"Einheitlishe feldstheorie von Gravitation und Elektrizitat" , trying to obtain this unifying 
goal, was criticized by Pauli and others, and Einstein died without success in this area. 



In Wheeler's book (IWheeler 19901 ). it is asserted that "... mass, going round and round 
in a circle, must produce a new kind of force, gravitomagnetism. Oersted discovered the 
creation of magnetism from moving electric charges in 1820, but gravitomagnetism has not 
yet been discovered, although it is the target of two great research enterprises (Gravity 
Probe B and LAGEOS project)". 

The term "gravitomagnetism" coined by Wheeler, is sometimes c alled "Gravitoele c- 



tromagnetism" and abbreviated as GEM by some recent authors (see iMashhoon 20031 ). 
The recent trend is to refer to it only as a set of formal analogies between Maxwell's field 
equations and an approximation to the Einstein field equations for general relativity, valid 
only for slowly moving particles. 

In this work, GEM will be given a different interpretation, not depending upon 
Einstein's general relativity approach to gravitation, and where Wheeler's force, plus other 
forces postulated by us, will be allowed to act not only upon mass particles, but also upon 
light quanta and conversely. 



We shall motivate this approach with additional historical considerations. As early 
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as 1918, Schrodinger (jShrodinger 19181 ) made a bright objection to Einstein's General 
Relativity handling of the stress-energy tensor, to which Einstein finally answered as follows: 
"... while sympathizing with these concerns, I am nonetheless convinced that a definition 
of the energy compon ents of the grav itational field more appropriate than the one I have 



given is not possible" (lEinstein 19181 ). No further comments 



Again, even as late in Einstein's life as 1950, he was cautio us about his general theory 



of relativity and unified field theory, and writes (lEinsten 19501 ) "The editors of Scientific 
American have asked me to write about my recent work which has just been published. It is 
a mathematical investigation concerning the foundations of field physics and also "As 
for my latest theoretical work, I do not feel justified in giving a detailed account of it before 
a wide group of readers interested in science. That should be done only with theories which 
have been adequately confirmed by experience" . These are the words of a great scientific 
mind indeed, and set an example to be followed, notwithstanding the fact that a theoretical 
work may precede its experimental confirmation, as it happened with Dirac's theoretical 
prediction of a positron, among examples of this fact taken from "modern physics" as 
distinguished from "classical phy sics" in the sense given to these terms by L. De Broglie, 



the founder of Wave Mechanics (jPe Broglie 19391 ). Besides, De Broglie emphasizes that 



modern physics is based on experiment, "constituting a method which deliberately produces 
certain given con ditions in order to see what phenomena are caused by these conditions" 
Jpe Broglie 19391 ). 



More precisely, De Broglie (1939) sets up a four steps procedure to obtain a sound 
physical theory: 

i) The experimental method allows certain laws to be established. 

ii) Theory then interprets these laws by establishing a connection with a single principle. 

iii) Theory uses the principle thus formulated in order to predict quantitatively other 
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phenomena. 

iv) Finally, recourse is had once again to experiment to verify the precision of the predictions 
made by theory. 

We shall adhere to this guiding procedure in relation with our proposed model of vector 
gravitational waves in the final section of this article, and shall also refer to "experiments" 
already considered in this area. 

In his 1937 book, De Broglie, having begun with a Chapter entitled "A general survey 
of present-day physics", in Chapter IV, section 4 "Relativity and Quanta" writes: "In its 
earlier form, the Theory of Relativity was concerned with the space and time co-ordinates 
only of observers in an state of uniform motion in a straight line. Later it was generalized by 
Einstein himself, and has provided a des cription of Gravit ation. This generalized Relativity 



however, I shall here leave on one side" (jPe Broglie 19391 ). 



In another book, published the same year in French (iDe Broglie 19521 ) De Broglie 
writes the following words to which we adhere: "But if the Special Theory of Relativity 
seems to be sufficiently confirmed by experiment, we believe it convenient to be a little 
less definite in that which refers to the Generalized Theory. The new phenomena which it 
predicts are in fact very small, and even though they may be observed, one can always ask 
whether they are really originated by the cause which Einstein's Theory attributes them or 
by some other perturbing cause, very small and negligeable while reasoning. Neither the 
very weak precession of Mercury's perihelion nor the very weak light rays deviation passing 
near the sun disk seem to supply an irrefutable proof of the exactness of the relativistic 
conception about gravitation: these phenomena do exist and have the order of magnitude 
predicted by Einstein's Theory, but his interpretation does not impose at all upon us." 

Besides Einstein's cautious words about his Generalized Theory already quoted, we 
might also mention the following episode concerning Feynman's first professional talk in 
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February 1941 held at the Princeton Physics Colloquium. Pauli objects to Feynman's 
lecture, and then Pauli said "Don't you agree, Professor Einstein?" , and Einstein's reply to 
Pauli was "no, the theory seemed possible, perhaps there was a conflict with the theory 
of gravitation , but after all the theory of gravitation was not so well established 
f lGleick imm . 

Recently, (July 2004) Stephen Hawking gave a talk on his new theory of quantum 
gravity at his department in Cambridge University. He used a complex mathematical 
technique called the "Euclidean path integral," first used by Richard Feynman, it has 
generally been applied to subatomic particles, but Hawking applies the idea to black holes 
(IPeplow 2004 ). 

Having thus reviewed what have been the main proposals for a theory of Gravitation, 
together with the mathematical and physical contexts which likely played a part in their 
formulation, we shall proceed to detail our proposal. 



2. Vibrations and Generalized Vector Waves 

The starting point h ere is the following obvious general consideration made by De 
Broglie (jPe Broglie 19391 ): "whatever magnitude it is that is traveling in the form of waves 
must pass through a periodic variation, the period itself being obviously equal to the time 
elapsing between the passing of two consecutive crests." These periodic variations result 
in a vibration at any given point. Thus we can define a physical wave as a vibration of 
some physical magnitude which travels through space-time with a definite velocity. Wave 
velocities are usually associated with characteristic values for hyperbolic systems of PDE's. 
Waves in continuous media whose origin is due to an initial disturbance imparted to the 
propagating media, and not to a vibration forced at some points of the media, also fit with 
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De Broglie's general definition of waves, since the initial disturbance can be interpreted as an 
initial energy provided to a continuously distributed system, formed by an infinite number 
of simple oscillators arranged as contiguous differential elements capable of vibration at each 
point of the media. Besides physical waves, we have to consider Schrodinger's "probability 
waves", and their origin is altogether different, so that the following considerations do not 
apply to them. Let us recall two of the most familiar and simple free vibrating systems, in 
order to define a generalized simple vibration system. 

a) The system mass-spring. If "x" denotes the horizontal displacement at time "t" of a mass 
"m" attached to a spring with stiffness "k", which moves without friction, and defining 
V = dx/dt, Newton's third Law of motion yields: 



m— + kx = (J 
dt 



Multiplication by "v" followed by time integration yields: 



—v'^ + —x^ = constant (2) 

b) A closed circuit formed by a coil or solenoid with inductance L and a condenser with 
capacitance C. Defining "g" as the electrical charge at a time "f ' stored by the condenser, 
and letting i — dq/dt stand for the electical current density, circulating in a coil with 
negligeable ohmic resistance, Kirchhoff 's Law yields: 



^di 1 



Multiplication by i followed by time integration yields: 
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-i^ + — — -q'^ = constant (4) 

We see, generalizing from these two classical examples, that the simplest physical 
systems (characterized by "lumped" parameters), in order to vibrate, require that they 
are able to store energy in two different forms [kinetic energy of the mass "m" and elastic 
potential energy of the spring for example a), magnetic energy associated to the magnetic 
flux of the coil, and electric energy associated to the condenser's electric field for b)]. 

So, if we have two scalar magnitudes Ui and U2, capable of representing two 
different energies when squared, and further related by U2 = dui/dt, they will vibrate 
when the physical law determined by the nature of the system is of the general form 
CiUi + C2{du2/dt) — 0, where Ci and C2 are two positive constants describing some lumped 
parameters, and this second differential equation when multiplied by dui/dt and integrated 
exhibits the "Conservation of Energy" Law for this system: Ci(tti)^ + C2('U2)^ = constant. 

De Broglie (1939) writes: "The importance of the part that certain physical magnitudes 
play result from the fact that there exists conservation theorems for these magnitudes". 
This is why the fundamental ordinary differential equation for our generalized system of 
two first order equations must be of the form —du2/dt proportional to Ui after defining 
U2 = dui/dt. Since ^(mi)^ = 2mi^(mi) = 2uiU2, this way of relating U2 and Ui ensures 
that when U2 vanishes, then the absolute value of Ui reaches a maximum in time, that is 
to say, at a time when one form of the system's total energy vanishes, the other form of 
energy reaches its maximum value, so that a total transfer has been made from one form of 
energy to the other, setting the basis for the oscillatory process. This basis is completed by 
the second relation between one magnitude and the time derivative of the other, namely, 
ui — —{Ci/ C2){du2/ dt) , which in view of the first one is equivalent to the conservation law: 
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— (Ml) + —[U2) = constant (5) 

The real positive valued constant of integration appearing in the right hand side, is 
the total energy of the system, which is acquired at some initial time to by imparting to 
the system initial values Mi(to) and ^2(^0) ioT the two scalar magnitudes Ui and U2- This 
in turn is accomplished through the action of some "forcing" terms, manifesting physical 
"forces" acting externally to the system considered during some time interval ending at to. 
If "b" denotes the horizontal force acting per unit mass ( "body force" ) along the positive 
X-axis in our mechanical example, and if "e" denotes the "electromotive force" set up in the 
coil by some externally induced field, then the differential equations for a) and b) become 
respectively: 



forced a) m— ^ + kx = mb; v = (6) 
at at 



forced b) L^ + (1/C)g = e; t = ^ (7) 

If the body force acts impulsively, with a large intensity "I" through a very short time 
interval centered at to, it can be mathematically idealized by means of the distribution 
16{t — to). The mass' horizontal displacement x{t) is continuous through the process, 
otherwise the spring and mass would break apart, but the velocity v{t), would appear as 
having a jump discontinuity at t = 0. 

Taking to = 0, and recalling that a piecewise continuously differentiable function, 
such as v{t), having a jump discontinuity at time t = 0, is a regular distribution, 
having a distributional d erivative composed of a regular part plus a singular part (see 



Miranda fc Pereyra 20071 ). but still satisfying the distributional version of the Fundamental 
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Theorem of classical Calculus: 

h 

J v'{t) dt = v{h) - v{-h) h positive , (8) 

-h 

we can compute the distributional definite integral from (-h) to h of both sides of the 
distributional equation defined by: 



— + [k/m)x = Id — > 
dt 

h h 

v{h) - v{-h) + {k/m) J x{t) dt = I , since j d{t) dt = 1 . (9) 

-h -h 

In the real, physical impulsive process, "/i" (defined so that v{t) = for t < —h) 
must be taken as a finite, though small, positive quantity In the idealized impulsive 
process, we take hmits as h tends to zero, thus getting v{0'^) — I as expected, as well as 
a;(0+) = x(0") = in this process. 

Another way of sustaining vibrations in such simple systems is by means of a forcing 
vibrating term (forced vibrations) of the form b{t) — M exp{iujt) , leading to resonance when 
cu^ = {k/m), "M" being the complex amplitude of the externally applied body force with 
angular frequency "a;" . 

These rather lengthy and known considerations about simple lumped parameters 
physical systems, whose free vibrations are governed by a system of two ordinary differential 
equations (ODE) of the first order of a very specific form tied to Conservations Laws, 
are justified as a natural way, leading from a discrete physical system to a continuous 
vector field distributed in three-dimensional space E3, by means of the analogy that can be 
established between the energy-densities of two vector fields within a small volume element 
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at a point of E3 and time "t", and the pair of magnitudes which vibrate in a lumped 
parameter system placed there at that time. 

Before going into the most general system, let us go back to our electrical circuit 

— * 

example. We know by Gauss' Theorem, that the electric displacement vector D ai a. 
surface point of the condenser's plate with a charge "g" uniformly distributed over the 
plate's surface with area A is such that D'^ — q^/A^, and this means that the vector D, 

— * 

and therefore the electric field E, normal to the plate surface, is proportional to "q" at any 
given instant of time "t". Now, since q = —LC{di/dt) and the mean magnetic induction B 
established by the current intensity "i" circulating through the coil is such that B — {L/A)i, 
we see that in the localized small region of space about a point x — {xi,X2,X3) where we 

— * 

assume our analog lumped parameter system to be concentrated, the electric field E is 
proportional to {—dB/dt) since q = —AC{dB/dt) = DA, and in the simplest case of a 
condenser formed by two plane parallel plates of area A and separation "W" , D/C — EW/A 
and -dB/dt = {W/A)E. Note that (W/A) has the dimension of 1/ length. 

In the same way, the other equation of the first order system, which is the definition 
i = dq/dt, together with B = (L/A)i and q ^ DA tell us that {l/L)B = dD/dt, 
and since we know that for a solenoid (B/L) has the dimensions of H/ length, we can 
write that the second, defining equation of the system can be written symbolically as 

— * — * — * — * 

dD/dt — {l/length)H, and the previous one as —dB/dt — {l/length)E. 

When we consider a distributed system, i.e., the fields E and H varying in space as 

— * — * 

well as time, it is to be expected that the ordinary time derivatives of D and B become 
first partial derivatives with respect to time, and that the right hand sides become some 
first order partial differential operator of H and E evaluated at time "t" as functions of the 

— * — * 

space variables x\, X2, xs, since only the net contributions of H and E are to be accounted 
for when they vary in small volume elements about point x — {xi, X2, X3). 
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The first order "d" operator appearing in the generahzed Stokes' Theorem for 
differential forms over smooth manifolds, is represented by means of the divergence, rot 
(curl) or gradient operators in particular cases. Now, the divergence operator takes vectors 
into scalars and the gradient takes scalars fields into vector valued fields. Thus, when only 
vector valued fields can be the domain and range of a first order "d" operator, the only 
remaining one is the rot operator, and we would naturally write as a system of first order 
PDE (that extends the above ODE first order homogeneous system) the following vector 
electromagnetic (VEM) equations: 



f = Vxi? (10) 



The above heuristic considerations are not meant to "prove" anything since the 
second equation of our PDE system is simply the differential form of Faraday's law: 
e.m.f. = —-^{magnetic flux). They are presented as a motivation for considering a 
Maxwellian type of system as a natural vector wave prototype or canonical form, at least 
when only two solenoidal vibrating vector fields are capable of storing the energy of the 
system. 

On the other hand, the distributed system stemming from the mass-spring lumped 
mechanical vibrating system, extends to the familiar distributed system known as "the 
vibrating string or membrane" , with a body force, such as gravity, now incorporated as a 
forcing term: — TAu + pb, with b — —gk, u{x, t) — u{x, t)k, x — xii, or x — xii + X2j, 
and A is the one or two dimensional scalar Laplacian, according to whether we are dealing 
with a string or a membrane. This is certainly a non homogeneous scalar wave equation 
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for the scalar unknown u{x,t), but in contrast with the previous electromagnetic system, 
now appear linear homogeneous boundary conditions at the finite boundary of the tensed 
string or membrane, becoming the first and simplest Quantum Model in physics ("T" 
stands for the horizontal initial tension of the string or membrane at rest, "p"is the 
mass per unit length or area at rest), first considered by Bernoulli and Euler yielding a 
discrete spectrum of allowable vibration frequencies associated to "stationary waves" . The 
associated eigenvalue problem and its solution was later generalized into the Sturm-Liouville 
Theory, as was anticipated by Poincare in the early 19*^^ hundreds. The higher dimensional 
generalization of the vibrating string equation is, for an isotropic homogeneous linearly 
elastic three dimensional body: 

P^ = V-T{u)+pb , (12) 

where T{u) is the stress tensor linearly associated with the strain tensor for the displacement 
u{x, t) at time "t" of a body point initially at x. The asymptotic distribution of the 
eigenfrequencies for a body of general shape is the subject of a classical paper by Weyl 
(1915). The elasto dynamics vector equation, 



/xAm + (A + /i)VV -M = p— -p6 (13) 



[where A and /i are Lame's constants, (IMiranda 19691 )]. admits two kinds of waves, which 



can be identified with the familiar P- Waves (compressional waves) and S- Waves (shear 
waves) of Seismology. P-Waves are not transversal, and this, together with its natural 
quantization under appropriate linear homogeneous boundary conditions, prescribed on a 
finite boundary manifold leading to Sturm-Liouville type of problems, make this distributed 
mechanical model substantially different from the distributed Maxwell model arising from 
the electric circuit model. This fact may lie at the root of the unsuccessful attempts to 



-17- 



explain light propagation by means of a hypothetical elastic carrier or "ether" , whose history 
has been extensively described by Whittaker (1973). Classical Maxwell-Lorentz VEM 
equations alone cannot explain quantized energy states. Besides, the elastic displacement 
vector field u{x, t) is NOT solenoidal, and neither the velocity field §iu{x, t) is solenoidal, so 
that Maxwell-Lorentz is not the appropriated canonical form for vector waves in this case. 

— * 

If we think of if = ■ui as a fundamental unknown vector field and the magnetic 

— * — * — * 

induction B — ^qH (in free space) as a secondary magnitude, and think oi E — U2 as 
fundamental vector field and the electric displacement D = eqE as secondary, we obtain the 
first order system: 



VxMi - eo-^ (14) 

V X 142 = ~^'^~df ■ 

In the same manner as we asked how a vibration could be initiated and sustained in 

— * — * 

time, it is natural to ask what can be the sources for fields U2 — E and Ui — H. 

In the case of the electric field, again Gauss' Theorem tells us that at points of space 
where an electric charge of volume density p is non zero, we must have V ■ U2 — p/ e^. If 
U2 — —VU with U — (l/47reor), then V • ■U2 = —AU — 5/eo, so that U2 is the electric 
field due to an electric monopole with unit charge at the origin of coordinates, since 
J J J S dxidx2dx3 — 1 for definite distributional volume integrals extended over closed 
domains having the origin as an interior point. 

So far, no magnetic monopoles have been found in nature, in spite of a prolonged 

— * 

search by Dirac and others, so that one must write always V ■ ui — everywhere. 
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A dditionally, due to the diff erential form of Ampere's Law for the "magnetomotive 
force" ( Bleaney fc Bleaney IQsJ ). moving charges constituting a vector density current j, 



act as an additional source term for V x besides, when we have to deal with conducting 
bodies, j = crE (Ohm's Law 



and when it is a matter with a "convection current" as 
Lorentz calls it (iLorentz 19091 ). j = pv, where v is the velocity of the freely moving charge 
passing through point x = (xi,X2,X3) at time "t". We shall not use Gaussian units as is 
usually done in most theoretical presentations, and this for several reasons, one of them 
being that when you do so, it veils the relation of the resulting wave speed to other 
important constants related to the propagating medium. 

We shall take eo = 8.8547 x 10-^^[C'^N-^m-^ and po = in x 10-'^[NA-'^], where 
we used the rationalized meter-kilogram-second-Coulomb system employed by Bleaney & 
Bleaney (1959). Coulomb's Law of force in free space has the form 

91(12^ , (16) 



47reo 

the force being repulsive when both charges have the same sign. 

When dealing with dielectric or magnetic media, two new constants are added, e 
and fi, which are multiplicative correcting factors, called the dielectric constant and the 
permeability of the medium respectively, so that for isotropic media, D = ee^E and 
B = hhqH. Of course, for free space is e = /i = 1. 

Finally, the corrected first order non- homogeneous system of PDE's for isotropic media 
in the presence of sources p and j adopt the following Maxwell-Lorentz Form (ML): 



Vxm = eeo^+j (17) 
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V X = ~^^o~df (18) 



V-ui = (19) 



V-iTa = — (20) 
eeo 



Writing the vector Laplacian A = VV • —V x Vx, we obtain as usual the following 
non homogeneous vector wave equations (VEMW): 



^""^^ = l^meo§^U2 + V (^) + /^/^o| (21) 

Aui = ///xoeeo— Ml - V X j (22) 
These two equations will be generalized to Vector Gravitational Waves (VGW). In free 

— * 

space j — p — and e — /j, — 1, and the speed of the vector waves turns out to be "c" , the 
speed of light. It is important to note that for us, "free space" is not necessarily identical 

— * 

to "vacuum" , and it will mean only a region of E3 where the source terms j and p vanish 
identically during the time interval of observation for the field disturbance propagation. 
Other than that, "free space" will not be further precised for the time being. 

The change in sign in the first order time derivatives of the first two ML equations is 
necessary for conservation reasons as we have pointed out. There is an asymmetry between 

— * — * 

the equations for ui and U2 with respect to the source terms: there is no "/' term for V x it 2 
in the second ML equation and there is no "'p' term for V • Mi in the third ML equation, 
but this asymmetry might not be intrinsic if magnetic monopoles should appear after all, 

— * 

and with them, a magnetic current j term in the second ML equation and this should be 
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kept in mind with possible generalizations of the ML system of equations to other pairs of 
vector fields. 

The ML system must be supplemented with the Lorentz force equation, expressing the 

— * 

force F acting on a charged particle with charge q moving with velocity v m a. region of 
space where both ui and U2 coexist: 



It should be noted that the source (or forcing) terms p and j in finite regions of space 
are not the only cause possible for setting the motion of vector waves according to equations 
VEMW, because one can think of sources acting "at infinity" . This rather vague way of 
speaking can be given a more precise meaning. 

It is known that analytic complex functions of a complex variable z can model 
bidimensional electrostatic fields. It is also known that negative powers in a Laurent 
series expansion about a singular point (which we shall take to be zero for simplicity) can 
represent the potential fields (equipotential curves and the orthogonal field fines of force) 
of the various multipoles through its real and imaginary parts, so that 1/z stands for an 
electric dipole oriented along the real axis, l/z"^ stands for an electric quadrupole, etc. But, 
what is relevant to our present observation is the following fact: the transformation of the 
complex z-plane into the complex w-plane defined hj w — 1/z carries negative powers of z 
into positive powers of w, and also when z tends to zero, the absolute value of w tends to 
infinity. 

If you look at the equipotential curves of the identity function f[z) = z, namely 
X — constant together with its field lines y — constant, with z — x + iy, x, y real variables, 
i=imaginary unit, then it is clear that such would be the field resulting from the physical 




(23) 
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field established between two parallel plane plates with charges equal in magnitude but 
opposite in sign, with the planes of the plates orthogonal to the x-axis, and getting farther 
and farther away to x = — oo and x = +00 respectively. This is what one would call "an 
electric dipole at infinity" without ambiguity whatsoever. Again, if you look at f{z) = z"^, 
and take two pairs of parallel plates with opposite electric charges, but one orthogonal to 
the y-axis, while the other pair is orthogonal to the x-axis, and let them recede to infinity 
in an obvious way, we would have obtained what one would call "an electric quadrupole at 
infinity", and so on for the higher order multipoles associated with the various powers of 
the variable "z". 

The gist of the above consideration is that, when we speak of a regular analytic 
function, such as f{z) = z'^, which we see free of singularities in the finite plane, in reality we 
are dealing with "a quadrupole at infinity" , so that regular analytic functions are deceptively 
"regular" in the sense that they hide their singular origin at infinity. Whenever dealing 
with the solution of field equations in unbounded regions, mathematicians have faced a non 



uniqueness situation concerning boundary va 



and Neumann problems of Potential Theory (IMiranda et al. 20071 ). and as physicists also 



ue problems such as the classical Dirichlet 



have to cope with time and again, coming up with some particular way of handling it, as 
it happens with Sommerfeld's Radiation Conditions for the Helmholtz equation in domains 
extending to infinity. In three space dimensions we have an analog of a Laurent series 
expansion with the Spherical Harmonics Expansions for solutions of Laplace's equation. 
Potential Theory is a common mathematical tool for both electrostatics and the classical 
Newtonian Theory of Gravitation, which, in the context of the Gravitomagnetic Theory 
which we shall discuss in this work, might aptly be called "Gravitostatics" , in spite of 
having been used to describe planetary motions with great success, which can be explained 
by the small error that ignoring gravitomagnetic effects would have introduced in those 
observed phenomena, but which cannot be ignored in other phenomena to be discussed 
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ahead. 

What has to be kept in mind from the above considerations is that "what happens 
at infinity" cannot be safely ignored when deahng with field equations in large regions of 
space. Besides, the idea of modeling cosmic rays as "singularities coming from infinity" 
cannot be discarded a priori. 



As a very simple and illuminating example taken from (iMiranda et al. 20071 ). of how to 
deal in a precise way with "multipoles at the infinity", consider Laplace's equation in polar 
coordinates to be solved under a Dirichlet boundary data "d" in the unbounded domain 
exterior to the unit disk centered at the origin: 



d ( dV\ d'^V 
r— r— + — - = for r > 1, < < 27r (24) 



dr \ dr 

The general solution, obtained by Fourier's method is 



O 

V{r, 0) = ^ log r + V {An cos n0 + 5„ sin n0) (C„r" + D„r-") + C (25) 



n=l 



with Q, C, An, Bn, Cn, Dn, arbitrary real constants. 

By taking Dn = —Cn we obtain an infinity of solutions V* vanishing for r = 1, which 
correspond to a source of charge Q located at the origin (or at the infinity) plus a series of 
multipoles of equal strength but opposite sign located at the origin and at infinity: 

V* (r, 0) = — log r +y {an cos ?7,0 + bn sin n0) (r" — r~") (26) 
2n ^ — ^ 

n=l 

Since = r"(cosn0 + i sinn0) and z""^ = r~"(cosn0 — i sinn0), after taking adequate 
real parts, the connection of the terms of our series with multipoles at the origin (represented 
by powers 2;"") and with multipoles at infinity (represented by powers z^) is apparent. 
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The conclusion is that in order to obtain uniqueness for the Dirichlet problem in 
the space of unbounded harmonic functions, it is necessary to prescribe the total electric 
flux Q for the field solution plus the complex amplitude (a^ — ihn) of the multipole z"' 
for each natural integer "n". Thus, in order to precise an unbounded solution, we must 
know beforehand what are the causes of unboundedness present at infinity in order that 
our final solution behaves asymptotically well for each natural "n". If the radiation from 
"outer space" comes only from a physical monopole but not from multipoles, then only 
Q must be prescribed in addition to the usual Dirichlet boundary data "(i", and we must 
take a„ = 6„ = for all natural "n". The final, unique solution is obtained in the form 
V = V* + Vd, where Vd is the harmonic function "regular at infinity" (i.e., bounded in two 
dimensions) such that 



Vd{l, (f)) = C + {An cos ncj) + Bn sin ncj)) = d{(f)) 



(27) 



n=l 



Asympt otically, V behaves like V* since V* is unbounded. Similar reasonings have been 



used in (jPower fc Miranda 19871 ) to obtain sufficient conditions for uniqueness in the case 



of exterior Stokes flows. 



Use was made in that paper, as well as in (jPower &: Miranda 19861 ). of the singular 
solution of Stokes' system of equations known as "rotlet", and this solution will be found 
useful for discussing angular momentum effects of our present The ory of vector gravitationa l 



waves. Full use was made of the "Oseenlet" and its multipoles in (IMiranda &: Power 19831 ). 
establishing existence and uniqueness in the case of the Dirichlet exterior boundary value 
problem for the homogeneous Oseen system, defining "Oseen hydrodynamic potentials" as 
analogues of classical single and double-layer surface potentials. 

Before ending our way of considering the establishing of classical Vector Electromagnetic 
Waves (VEMW) starting from necessary conditions to set up vibrations in a general way. 
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we may recall that the magnetic field H in fi:ee space generated by a constant current 
distribution that extends to infinity along a straight filament, with filamentary current /, is 
such that by Biot-Savart's Law, it is tangential at each point of a circle of radius r, centered 
at the filament and in plane normal to it, with a magnitude given by if = I/2'Kr. Thus 
in an idealized experiment where a physical plane normal to a thin long wire conducting 
/ was to be covered by iron filings, one would then see the filings move to form circular 
patterns. This image is an analogue of the vision we have had in mind in order to set up 
our proposed equations, but with a semi-infinite light ray instead of an infinitely long wire. 

Rather than anticipating, as Wheeler did, that a mass going round and round in a 
circle, must produce a new kind of force, with which we have no quarrel, we have taken two 
additional steps in order to define these new gravitomagnetic forces. 

Firstly, we have considered a constant steady mass current distribution along an 
infinitely long straight filament producing circular lines of the new force on planes normal to 
it, generating angular momentum in an obvious way. This certainly would be accomplished 
if one thinks of the Generahzed Gravitational Field as having two "components" satisfying 
the adequate analogue of Maxwell-Lorentz (ML), which might be called "gravitomagnetic" 
(Generalized Gravitational Ui) and "gravitoelectric" (Generalized Gravitational U2), with 
a minus sign preceding the source terms stemming from the fact that two positive masses 
attract in gravity, while two positive charges repel in electromagnetism as Mashhoon 
adequately observes. 

Secondly, and here is where the name "Lorentz-Compton approach" comes in, we 
allow the gravitomagnetic component to be generated not only by a mass current density 
j, but also by Light itself in the form of a current of photons incoming, as a semi- infinite 
light ray, into a "Trio-Hole" to be defined in section 4. If we had a monochromatic ray of 
frequency v say, then each photon would have an energy hv {h denoting Planck's constant). 
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and exactly as in the Compton effect, it could behave as a particle of equivalent mass 
hvj(? and equivalent momentum {hv/c)u, where "m" is a unit vector in the direction of 
the light ray, or normal to the light wave front if preferred. Highly energetic light beams 
as generators of matter and of the characteristics of the Created Universe, were considered 
in the or al Jewish tradit ions kept by Tzfat (Safed) learned men, and named "kav lines" 



by them (ILatapiat 20071 ) . They also considered the "Big Crunch" , a phenomenon prior to 
Gamow's "Big Bang", under the name of "tzim-tzum". Thinkers of the past have often 
had intuitions about our physical world that surprise us once and again. For a more recent 
historical example we can consider Descartes' 1625 Theory of vortices, whereby he declared 
that matter should not be made only of atoms, but also of a much smaller particle, which 
he conceived "tourbillon" (vortex). 

The constants cq, e, /iq, /i appearing in (ML) and (VEMW) will have gravitational 
analogs eoc,, e^, /iog; l^g- In free space we take eg = fig = 1, and eog will be related 
to Newton's constant G = 6.67259 x lO^^^m^ / Kg ■ sec^ by the obvious relation 
eog = l/47rG = 1.19 x lO^Kg ■ sec^/m^. 

Concerning fiQg, two options appear, one is to assume the speed of Vector Gravitational 
Waves described (VGW) below to be the speed of light "c" , and the other option would be 
to measure /iog directly and then compute the resulting wave speed as {1/ fiogeogY^'^. 

Adopting the first option with = 8.987044666 x lO^^rv? /sec^, we get /^og = 
9.3301241 X 10-'^'^ [m/ Kg]. 

With the previously assigned meanings for j and p, the proposed system (GEM) for 
the gravitomagnetic field Ui and gravitoelectric field U2 is: 



du2 , s 

VxMi = egeog—-j (28) 
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-i ^ dui 

V X M2 = -PgPOg-g^ (29) 



V-?2i = (30) 



V-U2 = — (31) 



And the Vector Gravitational Waves equations (VGW) are: 



AU2 = PgPOg^g^Og-Q^ y-^J - PgPOg-Q^ (32) 

Ami = figfiogegtog-^ + V X j (33) 

where j can be a photonic current vector density and not only a mass current density as in 
Mashhoon. 

The GEM system is not to be considered an approximation to geometric field equations 
like Einstein's, as Mashhoon views them, but on the contrary, GEM is regarded as an 
exact system replacing wholly Einstein's original geometric model concerning gravity. 
Gravitational waves appear as vector wa ves, solutions of exact linear wave equations, unlike 



Dirac's gravitational waves (jPirac 1960l ). which appear as approximations in the form of a 
D'Alembertian operator applied to each small magnitude h^i, resulting from the deviation 
of the fundamental tensor (7^1, with respect to the metric tensor associated to a totally 
fiat space-time when the Gravitation Constant is taken to be so small that the contracted 
Curvature Tensor R^^ can be approximated to zero. After neglecting terms of higher order 
in 7, Dirac obtains from Einstein's fundamental field equations the following: 
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dV dV 1 d"^ 



Dirac further chooses a Coordinate System so that the "harmonic conditions" = 
are satisfied, so that an inhomogeneous wave equation, expressed in this particular 
Coordinate System, is satisfied for each of the h^^s: 

Uh^^ = Wnjpf,^, with p^^ = T^^ - g^,uTx/2 (35) 

Finally, in domains where no matter is found, so that the material energy tensor T^j^ 
vanishes, Dirac gets Dh^u = as the equation to describe waves traveling in free space with 
the speed of Light "c" . 

In this paper, Dirac coins the word "gravitons" for particles of Gravitation's energy 
traveling with the speed of Light "c" , and states that the transfer from these Gravitational 
waves to Quantum Theory with the accompanying Hamiltonian Formalism is wholly not 
complicated ('ganz unkomplizie rt'). For details concerning nomenclature, we refer to Dirac's 



standard textbook (IDirac 19751 ). where it is explicitly recognized that ''f^ (the energy and 
momentum of the gravitational field) cannot be a tensor" (it is only a pseudo-tensor) and 
that "it is not possible to obtain an expression for the energy of the gravitational field 
satisfying both the conditions: i)when added to other forms of energy the total energy is 
conserved, and ii)the energy within a definite (three dimensional) region at a certain time 
is independent of the coordinate system. Thus, in general, gravitational energy cannot be 
localized. The best we can do is to use the pseudo-tensor which satisfies condition i) but 
not ii). It gives us approximate informatio n about grav itational energy, which in some 



special cases can be accurate". Dirac dixit! (IDirac 19751 ). Referring to the energy of Dirac's 



gravitational waves, Dirac states that "owing to the pseudo-tensor not being a real tensor, 
we do not get, in general, a clear result independent of the coordinate system. But there is 
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one special case in which we do get a clear result; namely, when the waves are all moving 



in the same direction." (iDirac 19751 ). In the face of all this, it is hard for us to see the 
need to cling to Einstein's General Theory of Relativity with its unnecessary complicated 
non-linear field equations in order to explain gravitational phenomena which can also be 
explained from a linear system of vector equations such as (GEM), which in addition 
predicts other phenomena in the astronomical scale. Even such things as the hypothetical 
"black holes" could be considered now within the framework of (GEM) with an analogue of 
strong laser beams (used for the source term j) employed by the 1997 Physics Nobel Prize 
winners Steven Chu, Claude Cohen- Tannoudji and William D. Phillips, in an experiment 
aptly described as "chilling with the light" , whereby we see a piece of matter contracted 
almost to a point and chilled near the absolute zero temperature. In our vision of the 
interaction of light and matter, this would be an earth lab scale analog of what could 
have been "the big crunch" , that is to say, such an initial singularity might not be due to 
gravitational collapse ensuing the end of the Universe expansion, but it could be due to 
light pressure concentrated at a point in E3 due to highly energetic incoming light rays 
"from infinity". Part of these light rays' energy would be reflected as outgoing light rays 
in a discrete number of directions, other part would be reflected as proto-matter, which 
besides traveling away from "the explosion center point", would be propelled to rotate 
about each ray alongside due to the gravitomagnetic field of the leading outgoing ray, which 
would then be transferring its initial energy after reflection to the rotating proto-nebulae as 
they move away, and such nebulae would be kept approximately in planes normal to the 
debilitating guiding reflected rays. We shall come back to this view, but in a mathematical 
way in section 4. 

We should not end our historical review without setting the record straight, giving due 
credit to H. A. Lorentz, who as far back as 1900, that is to say far ahead of Wheeler and 
Mashhoon, proposed the following, in his (originally written in Dutch) paper, which later 
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appeared in French ("Considerations sur la pesanteur" (iLorentz 19371 ). and from which we 
have freely translated to English, also using modern notation for the electric and magnetic 
fields and vector product, keeping our rationalized system of units, 

"It can be assumed that in the field of gravity there are two vectors E and B 
determined by equations (I), p being the density of weighable matter, and the force per unit 
mass shall be given by —r] ^^E + v x B^, r] being a certain positive coefficient. 

In every theory of gravitation it is important to examine which is the influence of the 
motion of heavenly bodies upon their mutual action. The answer to this question may be 
deduced from the preceding equations; besides, the problem is analogue in all points to a 
corresponding problem relative to the electromagnetic actions between charged particles. 
Etc., etc." 

Equations (I) referred to are given in page 208 of Lorentz' work (1937) under the form 
(in modern notation): 



V-D = p 



V-B = 



Vxif = pv + - 



p being the mass density at point x and time t, (recall that in Lorentz' electron theory, each 
particle with velocity v with respect to an ether assumed at rest, is likened to a small sphere 
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with a continuous density distribution inside it and vanishing at its surface boundary) . 

From Ecclesiastes 1:9 "The thing that hath been, it is that which shall be, and that 
which is done, is that which shall be done: and there is no new thing under the sun" . 

Thus we see that Wheeler's intuition and Mashhoon gravitomagnetic equations were 
long anticipated by Lorentz' genius. Probably, he didn't pursue this approach to gravitation 
because of his personal acquaintance and friendship with Einstein, whose towering figure by 
1916 may have inhibited him from opposing Einstein's views. 

But more important to us is that Lorentz' theoretical formulation of gravity as a strict 
analogue of electromagnetism, as it had been developed by 1900, in no way was tied to 
Einstein's theory of Gravitation, which came up in 1916. 

We follow his approach, without reference to an underlying ether of course, but besides 
a purely mass interaction, we have postulated a version of his system of equations in order 
to accommodate the force interaction between matter and light, analog to Compton's 
procedure, leading us a step beyond the simple mass-charge interchange with a sign change 
in the source terms. It is obvious that in Lorentz's formula the r) factor cannot multiply 
at the same time the gravito-electric and the gravito-magnetic force-field effects, since in 
the limit of the stationary case we would get Newton's force so that rj should be unity; so 
that this factor will be taken to multiply only the gravito-magnetic force term. We propose 
to obtain its value in such a way that the gravito-magnetic Larmor's precession formula 
yields the observed angular velocity of precession of Mercury's orbit. We have estimated 
the Sun's total angular momentum to be ^ 1.63 x 10^^ kgm^s^^ and have taken Mercury's 
orbital mean radius to be = 5.79 x 10^°m, to get uJuercury — 6.469 x 10~^^rad s~^ with a 
dimensionless 77 factor = 4.152 x 10^. 

A final word about Quantum Gravity is appropriate. Reading Hawking's speech at 
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the 17th International Conference on General Relativity and Gravitation, Dublin, given 
on July 21st 2004, it doesn't seem that his Theory could be considered as an established, 
definitive Gravitation Theory, pending on the Maldacena conjecture, and which further 
needs "regularizing" to cope with infinite action, as Hawking puts it. Thus, we feel that a 
theoretical space is opened for exact vector gravitational waves as solutions of a familiar 
linear system of PDEs, for which the modern quantum ele ctrodynamics procedure, already 
employed for Maxwell's equation (IFeymann &: Hibbs 19651 ). can now be considered for GEM 
equations, and this will be dealt with in Part II of the present work. Perhaps some might 
think that taking distance from Einstein's General Theory of Relativity would involve a 
denial of black holes. We think that it is necessary to distinguish between physical black 
holes, if there are such, and mathematical singularities (Schwartzchild black hole, Kerr 
black hole, etc,), which purportedly represent them, along with their observable properties. 
Different singularities may play the same role as those coming from Einstein's and related 
gravitational field equations. Our proposal includes singularities to be described in section 
4, and which we have baptized as "Trio-Holes" . 



3. The loss of energy in a binary star system due to gravitational waves. 

Since our vector gravitational wave equations are analogues to the corresponding 
electromagnetic wave equations, it is appropriate to recall some facts about electromagnetic 
radiation emission. In that case waves are formed by the time-change in the position and 
distribution of the electric charges in the system, and the same would be true when those 
"charges" are gravitational masses instead of electric charges. 

The simplest radiating systems are those represented by localized oscillating systems 
of charge and current density, as is done in classical electrodynami cs where the analysis 



is carried out on the basis of a multipole expansion (j Jackson 19751 ). Remembering that 
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the coefficients in that expansion are called multipole moments, and the total charge 
corresponds to the monopole moment, and since changing the total amount of charge in 
the system with time would violate the law of conservation of charge, then the way to 
accomplish wave-radiation is to have a time varying dipole moment, in the case of the 
electromagnetic wave. However in gravitational systems this is found to be impossible 
due to the violation of the law of conservation of angular momentum, so that one has 
to consider the next highest moment of mass distribution, the quadrupole moment, for 
possible emission of gravitational radiation, and this is coherent with the physics involved. 
The only known sources of gravitational waves strong enough to be detected by any known 
means are astrophysical phenomena, since only in deep space are massive bodies found 
moving fast enough to produce a sizeable wave signal. This phenomena includes binary star 
systems, pulsars, and many others. 

Binary star systems are, as their name implies, composed of two stars. The binary star 
system is one of the most common sources of gravitational radiation and it is the simplest 
case understood. Since our vector wave equations are formally the same as Maxwell's, 
the loss of orbital energy when the stars orbit each other about their center of gravity is 
explained mathematically in exactly the same way as is done in the similar case of orbiting 
electric charges, and this without invoking any General Relativity principles, or invoking 
any approximation of low velocities compared to the speed of light. This illustrates the 
relative simplicity of the proposed gravitational theory that conserves the principles of 
Special Relativity. 



4. The case of intense cosmogonic light beams and Trio-Holes 

We shall consider first an auxiliary quasi-steady photonic current density defined by 
j = dS in an otherwise free space, with 6 = S{r) = 6{x) ®6{y) ®6{z), where r = xi + yj + zk 
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as usual, a is a directional vector and (g) stands for the tensor product of one dimensional 
distributions, amounting to a "separation of variables procedure" for handling higher 
dimensional distributions when it is appropriate to do so. 

The associated auxiliary gravitomagnetic component ui would then, according to 
equation (VGW2), satisfy: Aui = V x j and V ■ u[ = 0, with j = a5. 

We are really interested in the gravitomagnetic field ul generated by a very narrow 
semi-infinite light ray, with a constant photonic current of intensity /, either outgoing from 
the origin to infinity along the positive z-axis direction, or incoming from infinity in the 
reverse direction towards the origin, where it is annihilated, which is the case we shall 
study below. The corresponding photonic current density j* can then be described in the 
following distributional form: 

/* = -Ikd{x) ® d{y) ® H{z) . (36) 

where H{z) is the Heaviside unit step function, with the known property H'{z) — S{z), so 
that j — —Ik5{x) (g) 5{y) (8) 5{z) — dj*/dz, and correspondingly, ui — dul/dz. 

Since H{z) — 1 for z > 0, J J j* ■ kdxdy — —I J J S{x) <S) S{y) dxdy — —I for 2; > 
when the domain of integration is any disk in a plane orthogonal to the cosmogonic beam 
and centered at the origin. However, we want to find u*i without recourse to a hypothetical 
analogue of Biot-Savart Law for differential elements of mass or photonic currents, which 
would be used, in order to know beforehand the tangential character of Ui by means of its 
circulation as in electromagnetism. 

For this purpose, we shall profit from what we already know about Stokes' fiows, and 
the "rotlet" singular solutions. The Stokes' system for a viscous incompressible fiow with 
unit kinematic viscosity and density is: 
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Av-^p=-b , (37) 

where b is the body force, and V • w = 0. . The fundamental solution of this system, 
associated to 6 = d6, has a pressure part, p = {a ■ r)/47rr^, and a velocity part S{d). Since 
V • V X t? = identically, it is readily found that, after computing the Stokeslet S{d) of 
vector strength a given by 



\- (a ■ r) — 



(38) 



[Fundamental solution of Ladyzhenskaya, see (ILadyzhenskaya 19641 )]. then the "rotlet" 
solution defined by rotlet{d) = V x (S'(a)), is such that ^{rotlet) = —V x (a5), and 
explicitly one gets from its definition that rotlet{d) = ^(a x r) (Power et al. 1987), and 
therefore we obtain that the auxiliary field ui = rotlet{d) = ^(a x r*) is the gravitomagnetic 
field generated by a photonic current density j = —d6. 

Finally, since Ui = du*i/dz, we can obtain u*i by means of an integration with respect 
to the z-variable after expressing r in cylindrical coordinates, f = pp + zk, and a = Ik, so 
that d X f= Ip6, P = p^ + z^. 



Observe that for z > 0, 



f*{x,y,z)= j j{x,y,t)dt = - j Ik5{x) ® 5{y) ® 5{t) dt = -Ik5{x) ® 5{y) ® H{z). (39) 

— oo — oo 

The expression for j*{x,y,z) is independent of "2;" for z > hj the above expression, 
but that is not the case for u*i{x, y, z) as we show below: 
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*/ N */ nN f9ul, f 2Ipd ^ 2ie 2/^' 

u^{x, y, z)-ui{x, y, 0) = y —dz = J ^ ^^^^^^ dz ^ — J {1+s ) ' ds 







p (l + s2)l/2 

(40) 



with s = z/p. The 9 direction explains the rotation of "proto-matter" (this term we 
shall leave undefined for the time being) about the 2;-axis. As z — > +oo with fixed p, 
u*i{x, y, z) —u*i{x, y, 0) = 219/ p asymptotically, as one would expect, since the semi-infinite 
photonic current tends to be seen, for large z > 0, as an infinite filamentary current with a 
1/p law for the associated gravitomagnetic field ui for all z. 

The solution u*i thus obtained is only approximate, being time independent, since due 
to the photonic annihilation assumed to take place at the origin, and depending upon the 
physical nature of the process taking place there (for instance, a total conversion of light 
into matter which would result at first from the simultaneous encounter at the origin of 
many incoming light rays), the mass density p at that point would be increasing with time, 
so that p{xi,X2,X3,t) = m{t)6{xi,X2,X3) and U2 will vary with time, with du2/dt different 
from zero, and now V x ui — —j + egeog{du2/dt) and Aui — —V x j + CgCogV x du2/dt. 
Since / // p{xi,X2.iX^,t) dxidx2dx^ — m{t) measures the total mass accumulated at 
the singularity (0, 0, 0) up to time "f ' , the magnitude of the error made by neglecting 
du2/dt will depend upon the time rate of change dm/dt of the increasing total mass m{t) 
accumulated at (0,0,0). 

Now V- {du2/dt) = —{1/egego) ■ [dp/dt] = — (l/e^ego) • {dm/dt) ■ 5{xi,X2,X3) which has 
an obvious solution given by 



du2 ^ ( 1 dm 1 \ 1 dm -i / 1 , , , 

• V I (41) 



dt \'^'KtgtQg dt r J Airegeog dt \r 



and 
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-L ( duo \ 1 dm -± ^ fl\ 

V X — ^ = • V X V - = 

\dt J A-Kegeog dt \r J 



—* — # —* — * 

so that it is still true that Aui — —V x j, but it will not be true that V x iTi = — j, since 



now 




(43) 



A better solution for the gravitomagnetic field (under the hypothetical conditions modeling 
our incoming photonic current) would be of the form {ul + u^*) , where 



This result for ul, even though approximate, obviously can be applied to explain the 
angular momentum that can be generated, by both incoming and reflected outgoing strong 
light rays ("kav" lines of Luria) along different a directions, in the cosmogonic processes 
after the Big Bang. The angular directions can be quantized in a more detailed Quantum 
Gravity model, and we shall present it soon in part II of the present work. 

It is also likely that this theory could have relevance in the atomic scale, related to 
some spectral lines features. 



In Einstein's geometrical approach to describe gravitational phenomena, it has been 
customary to baptize "Holes" according to mathematical singularities belonging to different 
metric tensors verifying the non linear field equations, and afterwards, the Cosmos has been 
searched for astronomical evidence of their "existence" . Our approach goes in the opposite 




(44) 



4.1. 



Trio-Holes 
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direction: we make hypotheses about possible physical processes that take place inside a 
given point, and then look after an adequate mathematical description of it, and then go to 
the Cosmos for corroboration. 

A singular point for our linear system of partial differential equations may be the siege 
for three different processes, namely: light beam annihilated and completely transformed 
into matter, light beam completely reflected, and light beam partially reflected and partially 
converted into matter. This is the reason to speak of "Trio- Holes" , and for simplicity, in 
this first paper, wc have worked out only the first process. Also, an exact time dependent 
solution for can be obtained by standard methods, but for qualitative purposes only, 
regarding the angular momentum present in the observed cosmos, our approximate solution 
is sufficient. 

Besides, the simplicity due to the linearity and familiarity of our proposed vector 
wave equations, allows the immediate treatment of complex situations by superposition of 
solutions. 

4.2. Universe mass previous to the Big-Bang 

We present a cosmological theory based in our view of gravitation as a form of 
interaction between matter and light, and in the equivalence of matter and energy, 
according to Einstein's equation E — mc^, taken in both directions, mass disintegration 
resulting into radiating energy, and conversely, matter resulting as a condensation of light, 
and also based in the behavior of light photons in collisions with matter as in the Compton 
effect. 

We postulate that "time" begins with a separation of a pervading hght, then occupying 
"infinity", from an "anti-light point" which becomes a "center". This "polar" situation 
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results into photonic beams incoming from "infinity" into the center. Initially, all of these 
photons are totally converted into mass at the center. As time goes by, mass increases, part 
of incoming light is reflected against the solid mass center, and part continues converting 
into mass in a "crunching" process, up to an instant when a critical magnitude is reached, 
when the big explosion known as "Big-Bang" occurs. Thus we can see matter in our present 
universe coming from different origins and ages, even before the "Big-Bang", as some 
astronomical observations and measurements presently suggest. 
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